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Abstract 

Most neutrino mass theories contain non-standard interactions (NSI) of neutrinos which can be 
either non-universal (NU) or flavor-changing (FC). We study the impact of such interactions on 
the determination of neutrino mixing parameters at a neutrino factory using the so-called "golden 
channels" ^v e — * for the measurement of #13. We show that a certain combination of FC 
interactions in neutrino source and earth matter can give exactly the same signal as oscillations 
arising due to #13. This implies that information about #13 can only be obtained if bounds on 
NSI are available. Taking into account the existing bounds on FC interactions, this leads to a 
drastic loss in sensitivity in #13, at least two orders of magnitude. A near detector at a neutrino 
factory offers the possibility to obtain stringent bounds on some NSI parameters. Such near site 
detector constitutes an essential ingredient of a neutrino factory and a necessary step towards the 
determination of #13 and subsequent study of leptonic CP violation. 
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I. INTRODUCTION 



From the long-standing solar |I| and atmospheric || neutrino anomalies we now have 
compelling evidence that an extension of the Standard Model (SM) of particle physics is 
necessary in the lepton sector. The simplest and most generic explanation of these exper- 
iments is provided by neutrino oscillations induced by neutrino masses and mixing. As is 
well known and accepted, the indication of the LSND experiment for oscillations at a 
large mass-squared difference can not be reconciled with solar and atmospheric data within 
a 3-neutrino framework. For recent four-neutrino Q analyses see Ref. || and references 
therein. For this reason we choose not to consider the LSND data and focus therefore on 
the simplest 3-neutrino scheme with the two mass-squared differences Am 2 ol < 10 -4 eV 2 
and Am 2 tm ~ 3 x 10~ 3 eV 2 ||. The corresponding lepton mixing matrix is parameterized 
by the three angles 9 12 , 023, #i 3 and one complex CP- violating phase 6 (see later for exact 
definitions) relevant in lepton-number-conserving neutrino oscillations 0, f§. It is known 
from atmospheric neutrino data that 623 has to be nearly maximal. On the other hand, 
data from present solar neutrino experiments favor a large value for the angle 9 12 P, 10|. 



An improved determination is expected from solar neutrino data and/or the results of the 
KamLAND experiment [ITT]. The value of the third angle #13 is not known at present, there 



is only the bound 

sin 2 2fl 13 <0.1 at 90% CL (1) 



implied by combining the results of the reactor experiments |L2| CHOOZ and Palo Verde with 
atmospheric data. Together with the requirement of solar neutrino oscillations this bound 
implies that sin 2 # 13 has to be small. The value of the phase 5 is completely unknown. 
Currently a new generation of long-baseline neutrino oscillation experiments using a neu- 



trino beam originating from the decay of muons in a storage ring is being discussed |T3 
These so-called neutrino factories are considered as the ideal tool to enhance our knowledge 
about neutrino mixing parameters. Besides the possibility to explore CP violation in the 
lepton sector an important aim of a neutrino factory will be a precise determination of #13; 
it is claimed that a measurement of #13 down to values of a fewxlO -4 will be possible fT4|| . 

In a large class of models beyond the Standard Model non-standard interactions (NSI) 
of neutrinos with matter arise. The simplest NSI do not require new interactions beyond 
those mediated by the Standard Model electroweak gauge bosons: it is simply nature of 
the leptonic charged and neutral current interactions which is non-standard due to the 
complexity of neutrino mixing 0. On the other hand NSI can also be mediated by the 
exchange of new particles with mass at the weak scale such as in some super-symmetric 
models with R-parity violating 16] interactions. Such non-standard flavor- violating 
physics can arise even in the absence of neutrino mass \ T7[ [18|j and can lead to non-universal 
(NU) or to flavor-changing (FC) neutrino interactions. 

Non-standard interactions of neutrinos affect their propagation in matter and the mag- 



2 



nitude of the effect depends on the interplay between conventional mass-induced neutrino 
oscillation features in matter [1E5 and those genuinely engendered by the NSI, which do 
not require neutrino mass [20|. Correspondingly their implications have been explored in 

1 El, atmospheric neutrinos 



a variety of contexts involving solar neutrinos [19, 20, 21 



24 , p5| , |26[| , other astrophysical sources |27], and the LSND experiment |29], The 



impact of non-standard interactions of neutrinos has also been considered from the point 
of view of future experiments involving solar neutrinos [31] as well as the upcoming neu- 
trino factories [32, 33]. Various aspects of NSI for a neutrino factory experiment have been 
considered in Refs. 0, g 0, [§, 0, |3j| |5| . 

In this paper we will consider from a phenomenological point of view the impact of NSI 
on the determination of neutrino mixing parameters at a neutrino factory. In particular 
we will focus on the { u e — > channels, which are supposed to be the "golden channels" 
for the measurement of #13. We extend our previous work |33| by taking simultaneously 
into account neutrino oscillations and the effect of NSI in neutrino source, propagation and 



detection [ 38] . We will show that a certain combination of FC interactions in source and 
propagation can give exactly the same signal as oscillations arising due to #13. This implies 
that information about #13 can only be obtained if bounds on NSI are available. In view of 
the existing bounds on FC interactions, this leads to a drastic loss in sensitivity in #13 , at 
least two orders of magnitude. 

All our considerations also apply to the determination of (9 13 in long baseline experiments 
using upgraded conventional beams . However, due to the different production pro- 

cesses involved, the NSI parameters relevant in the source of a conventional beam experiment 
differ from the ones at a neutrino factory. The same methods discussed here can be adapted 
to cover also that case. A detailed numerical consideration of conventional beam experiments 
goes beyond the scope of this work. 

The outline of the paper is as follows. In section [TT] we briefly sketch the theoretical 
motivation for NSI in the context of gauge theories of neutrino mass. In section |T| we discuss 
examples of low energy four-fermion Hamiltonians, which lead to NSI in neutrino source, 
propagation and detection. In section IV we review some bounds on NSI parameters obtained 
in the literature. In section [V] we present the framework of our numerical calculations 
and discuss the appearance rate in the presence of NSI and oscillations. In section [VT| 
we derive analytical expressions for this rate and formulate the oscillation-NSI confusion 
theorem. In section |V1I| we describe the simulation of a neutrino factory and our statistical 
method to investigate the possibilities of such an experiment to determine NSI and oscillation 
parameters. In section [VII]] we define sensitivity limits for sin 2 2# 13 and show our numerical 
results for the three baselines 700 km, 3000 km and 7000 km as a function of bounds on 
the relevant NSI parameters. We also discuss the sensitivity limits if information from two 
different baselines is combined. Finally we conclude in section [TX. 
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II. THEORETICAL MOTIVATION 



More often than not, models of neutrino mass are accompanied by NSI, leading generically 
to both oscillations and neutrino NSI in matter. The simplest are those NSI which arise 
from neutrino-mixing. The most straightforward example of this case is when neutrino 
masses follow from the admixture of isosinglet neutral heavy leptons as, for example, in 



seesaw schemes [^T|. These contain SU(2) <g> U(l) singlets with a gauge invariant Majorana 
mass term of the type Mr^vIv^ which breaks total lepton number symmetry, perhaps at a 
large 5*0(10) or left-right breaking scale. The masses of the light neutrinos are obtained by 
diagonalizing the mass matrix 

M r n 

(2) 



M L D 
D T M R 



in the basis z/, z/ c , where D is the standard SU(2) ® U(l) breaking Dirac mass term, and 
Mr = Mft is the large isosinglet Majorana mass and the M^vv term is an isotriplet 0. In 
left-right models the latter is generally suppressed as Ml qc 1/Mr. 

The structure of the associated effective weak V — A currents is rather complex J?|. 
The first point to notice is that the heavy isosinglets will mix with the ordinary isodoublet 
neutrinos in the charged current weak interaction. As a result, the mixing matrix describing 
the charged leptonic weak interaction is a rectangular matrix K @ which may be decomposed 
as 

K=(K L ,K H ) (3) 

where Kl and Kg are 3x3 matrices. The corresponding neutral weak interactions are 
described by a non-trivial matrix |7| 

P = K ] K . (4) 

In such models non-standard interactions of neutrinos with matter arise from the non- 
trivial structures of the charged and neutral weak currents. Note, however, that the smallness 
of neutrino mass, which follows due to the seesaw mechanism M vc q = Mr — DM^D 7 and 
the condition 

M L ^M R} (5) 

implies that the magnitude of neutrino NSI is expected to be negligible. However this 
need not be so in general. For example, since the number m of SU(2) <S> U(l) singlets is 
arbitrary, one may consider models with Majorana neutrinos based on any value of m. One 
can therefore extend the lepton sector of the SU{2) Cg> £7(1) theory by adding a set of two 



2-component isosinglet neutral fermions, denoted v c i and S{, in each generation [jT? |. In such 
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m = 6 models one can consider the 9x9 mass matrix 



PI 



D 
D T M 
M T jj 

[in the basis u, u c , S). The Majorana masses for the neutrinos are determined from 

M L = DM~ 1 nM T ~ 1 D T . 



(6) 



(7) 



In the limit fi — > the exact lepton number symmetry is recovered and neutrinos become 



massless |I7| . This provides an elegant way to generate neutrino masses without a super- 
heavy scale and automatically allows one to enhance the magnitude of neutrino NSI strengths 
by avoiding the constraints which arise from the smallness of neutrino masses presently 
indicated by the oscillation interpretation of solar and atmospheric neutrino data. 

The propagation of the light neutrinos is effectively described by a truncated mixing ma- 
trix Kl which is not unitary. This may lead to oscillation effects in matter, even if neutrinos 
were massless |P0fl . They maybe be resonant and therefore important in supernovae matter 
p0| , |27| . The strength of NSI is hence unrestricted by the magnitude of neutrino masses, only 
by universality limits, and may be large, at the few per cent level. The phenomenological 



implications of these models have been widely investigated [[43], [44 



An alternative way to induce neutrino NSI is in the context of the most general low- 
energy super-symmetry model, without R-parity conservation In addition to bilinear 
violation [[Tj], [|||] one may also have trilinear L violating couplings in the super-potential 



^ijkLiLjE k 

\jkLiQjD c k 



(8) 
(9) 



where L, Q, E c and D c are (chiral) super-fields which contain the usual lepton and quark 
SU{2) doublets and singlets, respectively, and i,j, k are generation indices. The couplings 
in Eq. ([8]) give rise at low energy to the following four-fermion effective Lagrangian for 
neutrinos interactions with <i-quark including 



£ e ff = ~2V2G F 22 VLal^VLp d R ^d R a,/3 = e,/i,T. 



(10) 



a,/3 



where the parameters £ a! 3 represent the strength of the effective interactions normalized to 
the Fermi constant Gp- One can identify explicitly, for example, the following non-standard 
flavor- conserving NSI couplings 



E 

3 

E 



I \' I 2 



4V2G F m 2 d 

I X' I 2 

^G F m\ L 



(11) 
(12) 
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and the FC coupling 



where rriq. L are the masses of the exchanged squarks and j = 1,2,3 denotes g?l, sl, & l, respec- 
tively. The existence of effective neutral current interactions contributing to the neutrino 
scattering off <i-quarks in matter, provides new flavor-conserving as well as flavor-changing 
terms for the matter potentials of neutrinos. Such NS1 are directly relevant for solar and 
atmospheric neutrino propagation Jj3 . 



Clearly, such neutrino NS1 are accompanied by non-zero neutrino masses. In fact one has 
a hybrid model for neutrino masses in which the atmospheric scale arises at the tree level 
while the solar neutrino scale is induced by loops which involve directly the NS1 coefficients 
in (HI). This way one obtains the co-existence of oscillations as well as NSI of neutrinos. The 
relative importance of NSI and oscillation features is parameter-dependent. 

An alternative variant of the above scheme is provided by some radiative models of 



neutrino masses such as the one discussed in |[49|| . In all such models NSI may arise from 
scalar interactions. 

Finally, we mention that unification provides an alternative and elegant way to in- 
duce neutrino NSI. For example unified super-symmetric models lead to NSI as a result 
of super-symmetric scalar lepton non-diagonal vertices induced by renormalization group 



evolution [Tj| |50 . In the special case of SU(5) the NSI may exist without neutrino mass. 
In 50(10) neutrino masses co-exist with neutrino NSI. 

In what follows we shall investigate the interplay of NSI-induced and neutrino-mass- 
induced (oscillation-induced) conversion of neutrinos at a neutrino factory and on how it 
can vitiate the, otherwise very precise, determination of neutrino oscillation parameters. 



III. EFFECTIVE FOUR-FERMION HAMILTON! ANS DESCRIBING NSI 

In this section we consider in some detail the simplest examples of effective low energy 
Hamiltonians, which lead to NSI in the source (S), propagation (P) and detection (D) of 
neutrinos in a neutrino factory experiment. 

In such an experiment neutrinos are produced by the decay of the stored muons [i + — > 
e + + v e + and the charge conjugated process. In the SM these processes are described by 
the effective four-fermion Hamiltonian 

H S SM = ^ [P„ (1 - 7 5 )7a /i] [e (1 - 7s)7 A *e] + h.c. (14) 

In addition to this SM term, resulting from the exchange of the W-boson, we consider now 
new processes fi + — > e + + v a + z/ M with any flavor a = e, fi, r for the neutrino related to the 
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positron. 1 We parametrize the corresponding NSI effective Hamiltonian by the coefficients 
e s : 

H s NSl = ^ [u, (1 - 75 )7a lA 4 a [e (1 - 7 5 )7 A v a ] + h-c. (15) 
Note that the value of the SM Fermi constant Gf is determined experimentally from the 



decay-width of the muon |51j, without measuring the flavor of the neutrinos. Therefore, we 
have the relation 

G^ = G F (\l + e s e f + J2 I4J 2 ) ' • (16) 

\ a=fj,,r / 

The term ef e leads to exactly the same final state as the SM process and must be added 
coherently, whereas ef and ef r lead to different final states and contribute incoherently to 
the decay. From Eq. ([16]) one learns that the high precision measurement of G G p P on its 
own (within an accuracy of 9 x 10~ 6 [J5TJ) does not constrain any of the parameters in the 
Hamiltonian Gp and ef a directly |5^| ; only the combination shown in Eq. (16) is constrained 
within the accuracy of the experimental measurement. 

The standard muon detectors under discussion for a neutrino factory experiment make 
use of charged current processes like + d — > ^ + u. The relevant effective Hamiltonian 
in the SM is given by 

H£ M = ( ^=[d(l- l5 ) lx u] [v, (1 - 75 ) 7 V] + h-c. (17) 

Here d (u) symbolizes any down- (up-) type quark. Similar to Eq. (|T5|) we consider the fol- 
lowing NSI four-fermion Hamiltonian: 



-| [d(l- l5 ) lx u] J2e^ [u a {l -75hV] +h.c. (18) 

The coefficients describe NU (a = jj) or FC (a = e, r) NSI in the detector, e.g. a non-zero 
leads to the process v T + d — > + u. 

In a long-baseline neutrino experiment a significant part of the neutrino path will cross 
the earth and hence neutrino NSI with earth matter must be taken into account. Let us 
consider the effective Hamiltonian describing the SM neutral current processes of neutrinos 
with a fermion / due to the exchange of the Z-boson v a + / — > v a + /: 



In this work we will consider the e — ► fi appearance channel and hence, we are interested only in the 
neutrino produced together with the positron (or the electron, for the charge conjugated processes). More 
generally, also processes /i + — » e + + v a + Dp with an arbitrary flavor combination (a, [3) = (e, /i, r) may 
occur. In this case the final state in the source can be different from the one in the SM. Such additional 
processes must be added incoherently to obtain the transition rate defined later in Eq. ( |26| ) . For simplicity 
we will not consider this possibility further. 
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n 



p _ Gp 

SM ~7I 



fiSv ~ ^7s)7a /] i 9 * t 1 ~ ^)7 A v a ] , (19) 



where g v and g^ A are the SM vector and axial couplings of the fermion /, see e.g. Ref. |)T| 
Sec. 10. In the SM this interaction is the same for all flavors and hence, has no effect on 
the propagation of the neutrino state - in contrast to the charged current interaction of 
u e with electrons. However, if NSI are present, we must also take into account processes 
v u + / — ► vp + / with arbitrary flavor combinations (aft) described by the Hamiltonian 

P _ G F 



f(g f v-g f A Khxf\ Mi -75)7%] • (20) 

a/3 



In Eqs. (15), (|18|), (ED) we have assumed for simplicity, that the NSI have the same V — A 
Lorentz structure as the SM interactions. This needs not to be the case in the most general 
extension of the SM involving, say, left-right symmetry, where many new NSI parameters 
can appear (see e.g. Ref. [J57L PSJ). However, the effects of NSI with V + A Lorentz structure 
are strongly suppressed since the left-right breaking scale should be rather high in order 
to account for the smallness of the neutrino masses indicated by solar and atmospheric 
experiments. Moreover, one expects that only certain combinations of parameters will be 
relevant for the experimental configuration we are considering here, and for any given theory 
our results can be mapped to the corresponding combination of parameters. 

Although different processes are relevant for source, propagation and detection, in a 
given model relations between the coefficients e^a may exist. However, such relations highly 
depend on the underlying model. In order to be model-independent we will treat all e^p as 
independent parameters. 



IV. BOUNDS ON NSI PARAMETERS 

In this section we review existing bounds on neutrino NSI obtained in the literature. The 
most direct upper bounds on the strength of NSI interactions arise from negative searches 
for neutrino oscillations pTI, [53|, |54| . For example the bounds on the transition probabilities 



Pu^vr < 3 - 4 x 1°~ 4 and Pv e ^v T < 2.6 x 10~ 2 obtained by CHORUS || yield the bounds 

| e CHORUS| < 1 g x 1Q -2 j | e CHORUS| < Q.16 . (21) 

With the superscript we indicate that the constrained quantity actually is a certain com- 
bination of NSI coefficients relevant in the neutrino source and detection processes for a 
given experiment (see later Eq. fl53p), which in general is different from the NSI coefficients 
relevant for neutrino factory experiments. 
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Recently, in Ref. |25| the strong evidence for oscillations of atmospheric neutrinos has 
been used to set upper bounds on NSI of neutrinos with the down quarks in earth matter: 



l< T l<3xlO- 2 , |ef r |<6x!0- 2 . (22) 



Similar bounds on the magnitude of neutrino NSI with electrons and up-type quarks may 



be derived 23 



Besides these direct bounds on neutrino NSI there is a lot of data constraining non- 
standard effects in charged lepton processes. However, it is very non-trivial in general to use 
these data to derive model-independent bounds on neutrino NSI coefficients. For a recent 
discussion see, for example Refs. |26], |29|, |35[. To obtain such bounds for neutrino interactions 
one must convert from the bounds on charged lepton processes making some assumption 



about SU(2)l symmetry. In Ref. [2^] the corresponding bounds for the charged leptons are 
multiplied by a factor of ~ 6.8, in order to take into account SU{2) L breaking effects. In 



this way the following bounds on neutrino NSI are derived from pure leptonic processes [22 : 



< 7 x 10 6 , |4|<3xl(T 2 . (23) 

As the neutrino production in a neutrino factory is also a pure leptonic process we take 
the bounds fl2"3] ) as order of magnitude estimates of the NSI at the neutrino source. From 
bounds on [i — > e conversion in muon scattering off nuclei and from those on flavour- violating 



hadronic tau decays the following bounds on neutrino NSI with quarks are derived p2 



I^|<7xl0- 5 , < 7 x 10- 2 . (24) 

We take this as an order of magnitude estimate for the NSI in propagation and detection at 
a neutrino factory experiment, since there also processes with quarks are involved. For the 



[a — t channel the bounds are of order [26 



< 5 x lO" 2 (25) 



and for the NU coefficients upper bounds of order 0.1 are derived. 

However, in Ref. []5S1 it has been stressed that in general no mo del- independent relation 
exists between NSI coefficients for charged leptons and neutrinos and only much weaker 
bounds of order 50% are derived using data from e + e~ colliders. This more conservative 
viewpoint has been exploited to show how FC neutrino interactions provide an excellent 
description of the solar neutrino data, while consistent with the oscillation description of 



the atmospheric data [23 
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V. THE APPEARANCE RATE IN A NEUTRINO FACTORY EXPERIMENT 



Let us consider the impact of NSI in source, propagation and detection on the e — > [i 
channel 2 at a neutrino factory. Starting from the decay of a fi + , we make the following 
ansatz for the rate at which a neutrino produced together with the positron leads to the 
production of a \i~ in the detector j38| : 

2 



n 



e/j, 



(26) 



V A s A p A 1 

a/3 

and similar for the charge conjugated process. Here we define the amplitudes describing the 
neutrino source and detection process as 



A x 



Sat 



z af3 



for X = S,D 



(27) 



and the amplitude AZ describes the propagation of the neutrino state from the production 
point to the detector. This amplitude is obtained from the solution of a Schrodinger equation 
with the Hamiltonian 
1 



2E, 



Udmg (0, Ami,, Att&J ^ + diag(V, 0, 0) + ^ V f , 

/ 



(28) 



which takes into account neutrino oscillations and SM interactions as well as NSI with the 
matter crossed by the neutrino beam. Here E v is the neutrino energy and V = V2G F N e 



is the matter potential due to the SM charged current interaction [fL9|l , where N e is the 
electron number density. The last term in Eq. (p8| ) describes the NSI with earth matter. 
The sum is over all fermions / present in matter, and Vf e a a is the coherent forward scattering 
amplitude of the process v a + f — > up + f, where Vf = ^/2GpNf, with the number density 
of the fermion / along the neutrino path given by Nf. We define an effective NSI coefficient 
for the propagation by normalizing all contributions to the down-quark potential Vd- 

V f 



€ aP — 



E 



v ^ 



(29) 



Adopting a basis where the charged lepton mass matrix is diagonal, the unitary matrix 
U in Eq. (p8| ) relates the neutrino fields in the basis where the neutrino mass matrix is 
diagonal to the neutrino fields in the basis where the interaction with the SM W-boson is 
diagonal [53|. We parameterize this matrix in the following way 0: 





( 1 





^ 




f 


C13 





Sl3 ^ 




( cia 




o\ 


u = u 23 u 13 u 12 = 





C23 











1 







-e- 45 si2 


Cl2 







v° 


-S 2 3 


C23 j 




\ 


-Sl3 





Cl3 J 




V 








(30) 



2 For the sake of clarity we keep the oscillation- inspired terminology "e — > [i appearance channel" ("//—> 
disappearance channel"). We actually mean the production of a wrong-sign (like-sign) muon in the 
detector, respectively. 
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where Sy = sinfljj and Cy = cos%. 



In this paper we consider the following simplified scenario. First, we take all e^ g real 



and we assume that they are the same for neutrinos and anti- neutrinos as in p5| . Second, 
Eqs.(]23|) and (pi]) suggest that constraints on FC interactions in the e— /x channel are about 3 
orders of magnitude stronger than in the other channels. This motivates the approximation 3 

4 WC J«° for X = S,P,D. (31) 

Third, we neglect the solar mass-squared difference, which implies also that the angle 12 
and the phase 6 disappear 0. Then we are left with the following neutrino propagation 
Hamiltonian 

H v = U 23 U 13 diag (0, 0, A) + diag( V, 0, 0) + Vr e p , (32) 

where we have defined A = Am1 tm /2E U and r = Vd/V = Nd/N e with r « 3 in earth 
matter. In the Hamiltonian (32) a sign change of A is equivalent to a sign change of V, 



which interchanges the evolution of neutrinos and anti-neutrinos. Therefore, it is sufficient 
to consider only the case A > 0, assuming that the neutrino factory is run in both polarities. 

A detector close to the front end of a neutrino factory - a so-called near detector - can 
be a very powerful tool to constrain NSI |3P|, |57| . Such a detector has to be situated 
at a short distance (a few 100 m) from the production region of the neutrinos, such that 
no oscillations with Am^ tm or Amg 0l can develop and matter effects are negligible. In our 
formalism this means that A^p = S a p, and the transition rate relevant for a near detector is 
simply given by 



V^ 5 A D , 



(33) 



It is clear that a near detector cannot provide any model- independent information on efg, 
and only a combination of ef^ and e^p is constrained. 

Some remarks are in order. Although the general motivation for NSI is that these ac- 
company models of neutrino mass generation, in our following phenomenological studies we 
will restrict our attention only to total lepton number conserving NSI. While this will suffice 
to make our point, it will on the other hand greatly simplify our analysis. This happens 
because in this particular case it is possible to distinguish the neutrino produced together 
with the electron from the one produced together with the muon, if the detector can de- 
termine the charge of the muon. We insist, however, that in a generic theory for NSI and 
neutrino masses also lepton number violating processes PO, p3l| are expected, due to the 



3 Note that in Ref. |35j much weaker sensitivities on FCI in the e — pb channel for the source at a neutrino 
factory are derived, of the order of 1CP 3 . A non-vanishing ef„ would imply yet an additional source for 
the confusion of oscillations and NSI. 
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Majorana nature of neutrinos (see Sec. H). Such effects would be an additional source for 
the confusion of NSI and oscillations. 

The off-diagonal elements e^JX = S, P, D) with a ^ (3 describe FC, whereas the diagonal 
elements with a = (3 lead to NU. In Eqs. ( p6[) and (|33|) we consider only processes with 
the same final states (in source and detector) as in the SM case. If additional processes 
are present, with different final states, the corresponding amplitudes have to be added 



incoherently to the rate [38 



We want to stress that our numerical results are not restricted only to the NSI resulting 
from the four-fermion operators discussed in Sec. [TTI[ The results apply to all kinds of non- 
standard physics in source, propagation and detection in a neutrino factory experiment, 
which can be parametrized like in Eqs. (p6|), (p7|), (|32|). In general the NSI parameter 
combinations involved in the quantum mechanical evolution of the neutrino system are 
model-dependent functions of the parameters appearing in the Lagrangian of a given theory, 
as discussed in Sec. |TI]. 



VI. THE OSCILLATION NSI CONFUSION THEOREM 

In this section we present analytical approximations for the transition rate and we show 
that within our simplified scenario NSI can lead to exactly the same signal at a neutrino 
factory as expected from genuine neutrino oscillations due to #13. Taking into account 
additional parameters like Am^ or CP violating phases in the lepton mixing matrix or in 
NSI can only bring more serious complications for the determination of 8 13 WA . 



For the understanding of the physics relevant for the numerical results which we will 
present in the following sections it is useful to derive an analytic expression for the appear- 
ance rate Eq. (p6|). To this aim we assume a constant matter potential V and consider the 



terms containing the NSI parameters e^g and S13 as a small perturbation of the Hamiltonian 
and calculate eigenvalues and eigenvectors of Eq. (|32l) up to first order in these small quan- 
tities. Then the appearance rate Eq. (^) is of second order in s 13 and e£j (X = S, P, D) and 
we make the interesting observation that only the three parameters (si 3 , ef T , e^ T ) appear. 

This is a special feature of the e — > u. channel under the approximation ([31]) and can 
be understood from Fig. [I] (a), where we show schematically the various contributions to 
this channel. The thick lines indicate the SM processes in source and detection and the 
dominating u. <-> r oscillations due to atmospheric oscillation parameters. The thin lines 
show the leading contributions of small quantities, which involve only the parameters S13, ef T 
and e^ T . With very thin lines we show some channels which involve more than one small 
quantity, and hence do not appear up to second order in the appearance rate: we consider 
only the detection of muons, therefore in leading order no effects of NSI in the detector 
show up because of Eq. ([[!]). Similarly no FC effects in the r <-» u. channel show up, since 
transitions from e to r flavor already involve a small quantity, either S13 or . We also 
note that no NU coefficient e;5 appears in leading order |34|, R5[, explaining the lack of 
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Figure 1: NSI contributions to the (a) appearance and (b) disappearance channel in the approx- 
imation Eq. The thick lines indicate the dominating processes, whereas the thin lines show 
the leading contributions of the small quantities S13 and e^. The processes shown with very thin 
lines are double suppressed in small quantities. 



sensitivity of neutrino factory experiments to NU parameters [[32 
Let us introduce the abbreviations 



(34) 



Then the expression for the appearance rate is a general quadratic form in the variables 
S13, es and e P : 



Tl eAL ss A s 2 13 + B s 13 e P + Cel + D e P e 5 + E el + F s 13 e s 



'en ~ ■"■ -3 13 T j-> 013^^ T u cp 

with the coefficients 

5 = 4 s 2 3 c 23 r 



-v 



n _ a „2 2 2 A 

— ^ *23 c 23 ' A-V 



A+y • 2 (A-V)L , ■ 2 VI • 2 AL 

^7 sm — 2 — ^ s IT — s IT 



— 4 s 2 3 c 23 r A _ y 



V ■ 2 (A-V)L . ■ 2 VL V ■ 2 AL 
sin 2 (A-V> _ ^2 ^ _ (! _ 2 J) sin 2 AL 



E = 4sl 3 4 3 sm 2 ^, 



F = 4s 2 23 C2 3 ^ v 



■ 2 (A-V)L ■ 2 VL 1 • 2 AL 

sm i — ^ sm ^ + sm ^ 



(35) 



(36) 
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where L is the distance between neutrino source and detector. The appearance rate for anti- 
neutrinos TZgfi is obtained by replacing V — > —V in Eq. These analytic expressions are 
in agreement with numerical calculations within a few % in the relevant parameter range. 
In general all coefficients (A, ... , F) are of the same order of magnitude and depend on 
neutrino energy (via A), on the baseline and on the sign of V (neutrinos or anti- neutrinos) 
in a nontrivial way. 

Performing a similar consideration for the \x disappearance channel one finds that the 
rate TZ^, defined in a way similar to Eq. (|26|) , contains terms of all powers of the small 
quantities. This is illustrated in Fig. [1] (b). The zeroth order contribution corresponds to 
<-> r oscillations with atmospheric neutrino oscillation parameters: in contrast to Fig. p] (a), 
in the case of Fig. [I] (b) there is a channel involving only thick lines. However, we find that up 
to second order in small quantities only the parameters (e^ T , e^, e^ T , , ef r , , e^ ) appear. 
The important observation is that none of the three parameters (sis, ef r , e^ T ) relevant for the 
e — ► \i channel appears. Therefore, no additional information on these parameters can be 
obtained by considering the disappearance channel. An analysis of this channel including 
all the parameters listed above goes beyond the scope of this paper and we will not consider 
it any further here. 

Let us compare the transition rate for pure oscillations (eg = ep = 0) 

K s ;(s 13 ) = As 2 13 (37) 
with the transition rate without oscillations (S13 = 0) but non-zero NSI coefficients 

n^\e s , e P ) = C 4 + D e P e s + E e% . (38) 
With the expressions for A,C,D,E given in Eq. (B9) it is easy to check that if the relation 

e s = re P (39) 
holds, oscillations are indistinguishable from NSI. More precisely, we obtain 

< SI (^p,ep)=7^ c ( Sl3 ) (40) 

with 

s?3 = r24 i±^H. (41) 

This means that for each value of S13 there is a pair of NSI parameters (es, ep) determined by 
Eqs. ([!9|) and ( flip which in our approximation leads to exactly the same signal as oscillations 
due to S13. This includes both energy and baseline dependence, for both neutrinos and anti- 
neutrinos. We call this the "oscillation-NSI confusion theorem" . 

Of course, relation (|39| ) represents a fine-tuning of the parameters ep and es- However, 
as long as this relation cannot be excluded one has to consider this possibility. Moreover, 
in a realistic experiment with finite errors and statistical uncertainties there will be a re- 
gion around the point in the (es,£p) plane corresponding to Eqs. ( [3~9|) and flHf) in which 
oscillations cannot be distinguished from NSI. 
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VII. THE SIMULATION OF A NEUTRINO FACTORY EXPERIMENT 



In our numerical calculations we assume a neutrino factory with an energy of 50 GeV 
for the stored muons and 2 x 10 20 useful muon decays of each polarity per year for a period 
of 5 years. We consider a magnetized iron calorimeter with a mass of 40 kt. The neutrino 
detection threshold is set to 4 GeV, the energy resolution of the detector is approximated by 
a Gaussian resolution function with AE U /E U = 10% and we use 20 bins in neutrino energy. 
We do not include any backgrounds, efficiencies and errors in the particle identification. The 
amplitude A^a describing the neutrino propagation is obtained by numerically solving the 
neutrino evolution equation with the Hamiltonian fl32|), using the average matter density 



along each baseline. In Ref. |58| it was shown that this is an excellent approximation as 
long as the baseline is shorter than approximately 10 000 km, i.e. as long as it does not cross 
the core. Then the transition rate Eq. (|26|) is folded with neutrino flux, cross section and 



energy resolution function to in order obtain the expected event rates in the detector. For 



further details see Refs. 14 



Our "observables" are the event rates for the appearance channel n % v (np) for neutrinos 
(anti- neutrinos) in each energy bin i. We fix the atmospheric oscillation parameters 4 at their 
best fit values given in Am 2 tm = 3 x 10~ 3 eV 2 and sin 2 26*23 = 1- Hence, at a given 

baseline, the event rates depend on the three parameters Si 3 , es and ep, where we have 
introduced the abbreviation = sin 2 20 13 . 

In order to evaluate the impact of e$ and ep on the capability of a neutrino factory to 
measure S\% we proceed as follows. To test a given point (Sf 3 , e° s , e° P ) in the parameter 
space we calculate the event rates (nl) = n^S^, e° s , e° P ) with x — v, v and construct a x 2 
appropriate for a Poisson distribution 5 : 



X 2 (Si3, es, ep; 5° 3 , e° s , e° F 



pj 



EE 



( 71 

<(Ws,ep)-(<)° + (<)°ln 



ni(S 13 ,e s ,e P ) 



(42) 



Thus we obtain an allowed region in the three dimensional space V spanned by (S13, es, ep) 
in the usual way. Note that the minimum of the \ 2 defined in Eq. ( ^2]) is zero and occurs 
at (5*13, e s , ep) = (Sl 3 , e s , e P ). Therefore, the allowed region at the CL a is given by the set 
of all points in V which fulfill 

X 2 (5 13 , e s , e P ; S° 13 , e% e° P ) < A X 2 a , (43) 

where A% 2 is determined by a ^-distribution with 3 degrees of freedom. For simplicity we 
consider only starting values with e° s = e P = 0, i.e. pure oscillations. 



4 They will be determined at the neutrino factory with high accuracy from the disappearance channel. 

5 We are dealing with a counting experiment with eventually very low counts. 
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Figure 2: Allowed regions at 90% CL (black solid lines) in the S13 = plane for three different 
starting values 5? 3 = (3.2 x 1(T 3 , 1(T 3 , 1(T 4 ) and e° s = e% = always. The baseline is 3000 km. 
The gray lines indicate points with the same event rates as the starting point (gray solid for 
neutrinos, gray dashed for anti- neutrinos) . 



VIII. SENSITIVITY LIMITS FOR sin 2 26>i 3 

In Fig. Q we show the allowed regions in the S13 = plane for a baseline of 3 000 km 
and three different starting values for 5° 3 . In gray we show the lines with the same total 
event rates (solid for neutrinos, dashed for anti-neutrinos) as in the starting point. The 
general shape of these lines can immediately be understood from Eq. fl35|). The small 
regions delimited by the dark solid lines arise from a global fit procedure including also the 
information on the spectra with the expected energy resolution as described above. One 
notices that these confidence regions follow closely the intersection of the lines of constant 
rates for neutrinos and anti-neutrinos. This means that most information is obtained from 



simultaneously taking into account neutrino and anti- neutrino rates [S3|. This follows from 
the fact that the allowed regions extend as long as the lines of constant neutrino and anti- 
neutrino rate are close to each other, i.e. both rates are similar to the ones in the test point. 
We conclude that it is important to run the neutrino factory in both polarities. On the 
other hand we learn that most of the information is contained in the total rates; the spectral 
information is not very important: our results are rather insensitive to variations of the 
number of energy bins and of the energy resolution assumed. 

Solutions in the Si 3 = plane of the type as shown in Fig. |] always exist, irrespective 
of the starting value S® 3 . This is a consequence of the confusion theorem we have presented 
in Sec. [V|. However, the magnitude of the required NSI parameters and ep strongly 
depends on the size of S® 3 as can be seen from Eq. ( filD or Fig. |[ Thus, it is only possible 
to derive a limit on S13 if there is a limit on es and/or ep; a neutrino factory can only test 
a certain value of S13 if the values of e^ and/or ep, which lead to the same signal, are ruled 
out by some other measurement. In Fig. ||| we show the attainable sensitivity limit on S13 
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Figure 3: Sensitivity limits at 90% CL on sin 2 2^13 attainable if a bound on ep (left panel), e 2 s 
(middle panel) or ep + e 2 s (right panel) is given. The dotted line is for a baseline of 700 km, the 
dash-dotted for 3 000 km and the dashed line for 7000 km. The horizontal black line illustrates 
the order of magnitude of current limits on the NSI parameter in order to "guide-the-eye" . The 
vertical gray band shows the range of possible sensitivities without NSI The diagonal solid 
line is the theoretical bound derived from our confusion theorem. 

as a function of different limits on the NSI parameters (the present estimated NSI limits are 
indicated by the solid horizontal lines). We define this sensitivity on S13 in the following 
way. For a dense grid of starting values S"° 3 in the range 10~ 6 — 10 _1 we calculate the 90% 
CL allowed region in the S\% = plane as shown in Fig. 0. For each value of S® 3 we show 
the minimum value of ep (left panel), e| (middle panel) or ep + e| (right panel) inside this 
allowed region. The neutrino factory is sensitive to this value of S13 only if there is a bound 
on this NSI parameter (combination of parameters), which is smaller than this minimum 
value. 

Any experiment (e.g. like a near detector at a neutrino factory) trying to measure es 
and ep will only restrict a certain combination of the NSI parameters in the source and the 
detector used in this particular experiment (see Eq. ([33])). In general it will be very difficult 
to translate such a result into a bound on es, and even more difficult on ep, in a model- 
independent way. It is however to be expected that the constrained combination depends on 
the square of the NSI parameters since the transition rate TZ oc e 2 . Therefore we show the 
results for the three simple functions e|, ep and e| + ep as mentioned above. The left hand 
panel of Fig. |3] shows the sensitivity limit if there is a bound on ep and all values of es are 
allowed. It seems very difficult to obtain such a bound in a model-independent way, because 
it is not possible to probe directly the NSI parameters relevant in neutrino propagation. As 
recently stressed in Ref. in general it is not possible to use bounds on similar processes 
involving charged leptons. Moreover, many different processes may contribute to ep (see 
Eq. (p9"D). In the middle panel we show the sensitivity limits for a bound on e| and all 
values allowed for ep, which is probably the most realistic case because it should be possible 
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Figure 4: Allowed regions at 90% CL (black solid lines) in the S13 = plane for three different 
baselines (L = 1000km, 3 000km, 5 000km) and S% 3 = 10~ 3 , e° s = e° p = 0. The gray lines indicate 
points with the same event rates as the starting point (gray solid for neutrinos, gray dashed for 
anti- neutrinos). The right hand panel is a blow up of the left hand panel. 

to constrain €3 with a near detector setup. In the right hand panel we display the optimal 
situation, if a bound on the combination e| + ep is available. 

The diagonal solid line in Fig. |3] shows the theoretical bound implied by the oscillation- 
NSI confusion theorem Eqs. (^) and (fH|) . This bound corresponds to the best possible 
situation, which can be achieved only for large values of S13 due to large event numbers. For 
smaller values of Si 3 the realistic bound gets worse because of statistical limitations due to 
small event numbers. The numerical differences between the three plots in Fig. |] is due to 
the different symmetries of the used function of e$ and ep with respect to the symmetry of 
the allowed regions in the S13 = plane. For small values of e| + ep (right panel) the bounds 
converge to the sensitivity limits obtained without taking into account NSI |I4[ . The range 
of these limits for the three different baselines is shown as the gray vertical band. 

We can understand the behavior of the sensitivity limits in Fig. |3] by considering the 
allowed regions in the S13 = plane for different baselines, as shown in Fig. [|. For small 
baselines the allowed region is roughly a circle. Therefore a bound on an individual ep or 
€3 is useless for a sensitive determination of Si 3 (see left and middle panel of Fig. |). We 
conclude that if only a baseline of 700 km is available it is mandatory to establish solid 
bounds on both NSI parameters. However, a bound on ep + e| is most suitable for small 
baselines and in this case L = 700 km can do even better than longer baselines (right panel 
of Fig. H). For longer baselines the allowed regions in the Si 3 = plane become smaller (see 
Fig. [D and hence, also a bound on an individual NSI parameter is useful. 

With the solid horizontal lines in Fig. |] we illustrate the order of magnitude of current 
bounds on NSI parameters. For ep we show the bound given in Eq. ([24]), while for the es 
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Figure 5: Sensitivity limits at 90% CL on sin 2 29\3 attainable if a bound on e 2 p (left panel), 



e| (middle panel) or 



+ ei 



(right panel) is given. The dotted line is for the baseline 



combination 700 km & 3 000 km, the dash-dotted for 700 km & 7 000 km and the dashed line for 
3 000 km & 7000 km. The horizontal black line illustrates the order of magnitude of existing limits 
on the NSI parameter. The vertical gray band shows the range of possible sensitivities without 



NSI [14|. The diagonal solid line is the theoretical bound derived from the confusion theorem. 



case we use the bound given in Eq. (|23|). It is clearly visible that using even these rather 
optimistic bounds on the relevant NSI parameters the sensitivity of a neutrino factory is at 
its best S13 ~ 1CT 3 - compared to S13 ~ 1CT 5 - in the absence of any NSI. The sensitivity 
is deteriorated by two orders of magnitude for all baselines. Let us stress again that we are 
not using the bounds from Eqs. fl24"D and (p3|) in our analysis because they are derived under 
some model-dependent assumptions from non-neutrino processes. The horizontal lines in 
Fig. |3] should merely "guide-the-eye" in reading the plots, they simply give a rough idea of 
the order of magnitude of existing bounds. 

The geometry of the currently discussed muon storage rings offers the striking possibility 
to illuminate two detectors at different baselines with neutrinos from one neutrino factory. 
With this in mind, let us investigate to which extent the sensitivities for S13 can be improved 
by combining the information of two baselines. From Fig. (| we find that although the shape 
of the allowed regions in the S13 = plane is very different for different baselines, they all 
have a common intersection. Moreover, considering the lines in the S13 = plane which 
have the same event rates as the starting value S® 3 (shown as gray solid lines for neutrinos 
and gray dashed lines for anti-neutrinos) we observe that all iso-rate lines meet in a single 
point. This is again due to our oscillation-NSI confusion theorem: Eqs. (|39|) and fl4l"|) do not 
depend on the baseline. Thus even the combination of baselines cannot lift this degeneracy. 

In Fig. [| we show sensitivity limits for S13 obtained as before, but now we use the sum of 
the x 2 -functions Eq. (H2T) for two different baselines. In comparison with Fig. R] we observe 



that the theoretical bound can now be achieved somewhat easier. However, in contrast 



to the case considered in our previous work [33 , where only NSI with the earth matter 
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are taken into account, in our present more realistic situation including also effects in the 
neutrino source even a combination of baselines does not resolve the confusion problem. 
With the current bounds on the relevant NSI parameters the sensitivity is Sis ^ 10 -3 , 
which coincides with the sensitivity obtained at a single baseline. Again this sensitivity is 
two orders of magnitude worse than without NSI. 



IX. CONCLUSIONS 

In this paper we have considered the impact of non-standard neutrino interactions on the 
determination of neutrino mixing parameters at a neutrino factory. In particular we have 
focused on the so-called "golden channels" for the measurement of #13, namely the v e — > 



channels. We have extended our previous work |33| by taking into account both the effects 
of neutrino oscillations as well as the effect of NSI at the neutrino source, propagation and 
detection ||38|| . Within a very good approximation we have explicitly demonstrated how a 
certain combination of FC interactions in source and propagation can produce exactly the 
same signal as would be expected from oscillations arising due to #13 . This implies that 
information about #13 can only be obtained if bounds on NSI parameters are available and 
that all one can achieve at a neutrino factory is a correlated oscillation-NSI study. In view 
of the current estimates of the bounds on FC interactions, this leads to a drastic loss in 
sensitivity in # 13 , at least two orders of magnitude. 

In order to improve the situation it is mandatory to obtain better bounds on e$ and 
ep at the e ~ 10 -4 — 10~ 3 level, which is several orders of magnitude more stringent than 
current limits. This unexpected complication should be taken into account in the design 
of a neutrino factory. On the other hand, a neutrino factory may also offer the possibility 
to obtain these very stringent limits on the NSI parameters. Using a small near detector 
(L ~ 100 m) with very good particle identification for taus it would be possible to restrict the 
Ti eT down to 10- 8 -10" 6 . Since TZ PT oc e 2 this translates into a bound for e ~ 10 4 — 10 3 |E7| 



Thus the near site physics program of a neutrino factory is a necessary and very important 
part of the long baseline program. 
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